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ABSTRACT 

Aims. The relative cosmic variance (cr„) is a fundamental source of uncertainty in pencil-heam surveys and, as a particular case of 
count-in-cell statistics, can be used to estimate the bias between galaxies and their underlying dark-matter distribution. Our goal is 
to test the significance of the clustering information encoded in the cr^ measured in the ALHAMBRA survey. 

Methods. We measure the cosmic variance of several galaxy populations selected with 6-hand luminosity at 0.35 < z < 1.05 as the 
intrinsic dispersion in the number density distribution derived from the 48 ALHAMBRA suhfields. We compare the ohservational (Ty 
with the cosmic variance of the dark matter expected from the theory, cr^ dm- This provides an estimation of the galaxy bias b. 
Results. The galaxy bias from the cosmic variance is in excellent agreement with the bias estimated by two-point correlation function 
analysis in ALHAMBRA. This holds for different redshift bins, for red and blue subsamples, and for several 6-band luminosity 
selections. We hnd that b increases with the 6-band luminosity and the redshift, as expected from previous work. Moreover, red 
galaxies have a larger bias than blue galaxies, with a relative bias of fcjei = 1.4 ± 0.2. 

Conclusions. Our results demonstrate that the cosmic variance measured in ALHAMBRA is due to the clustering of galaxies and 
can be used to characterise the cr^ affecting pencil-beam surveys. In addition, it can also be used to estimate the galaxy bias b from a 
method independent of correlation functions. 

Key words. Galaxies: statistics - Cosmology : dark matter 


1. Introduction 

One fundamental uncertainty in any observational measurement 
derived from galaxy surveys is the relative cosmic variance (cr^, 
also called sample variance), arising from the underlying large- 

* Based on observations collected at the German-Spanish Astronom¬ 
ical Center, Calar Alto, jointly operated by the Max-Planck-Institut fur 
Astronomie (MPIA) at Heidelberg and the Instituto de Astrofisica de 
Andalucia (CSIC) 


scale density fluctuations and leading to variances larger than 
those expected from simple Poisson statistics. The most efficient 
way to tackle the cosmic variance is to split the survey into sev¬ 
eral independent areas on the sky. This minimises the sampling 
problem and it is better than increasing the volume in a wide 
contiguous field (e.g. Driver & Robotham 2010). However, the 
uncertainties in many existing surveys are dominated by the cos¬ 
mic variance because of the observational constraints (depth vs 
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area). Thus, a proper estimation of cr^ is needed to fully describe 
the error budget in deep cosmological surveys. 

The impact of the cosmic variance in a given survey and red- 
shift range can be estimated using two basic methods: theoret¬ 
ically, by analysing cosmological simulations (e.g. Somerville 
et al. 2004; Trent! & Stiavelli 2008; Stringer et al. 2009; Moster 
et al. 2011), or empirically, by sampling a larger survey (e.g. 
Driver & Robotham 2010; Lopez-Sanjuan et al. 2014; Keenan 
et al. 2014). Unfortunately, the value of previous empirical work 
is limited by the precise understanding of the measured cosmic 
variance: the estimated cr^ must encode, by definition, the clus¬ 
tering of the studied populations, but no previous work in the 
literature has tested this requirement. 

The galaxy bias b is the relationship between the spatial dis¬ 
tribution of galaxies and the underlying dark-matter density field 
(Kaiser 1984; Bardeen et al. 1986; Mo & White 1996). The 
cosmic variance cr^ is a particular case of count-in-cell statis¬ 
tics (Peebles 1980; Efstathiou et al. 1990; Efstathiou 1995; An- 
dreani et al. 1994; Adelberger et al. 1998; Dekel & Lahav 1999; 
Marinoni et al. 2005; Robertson 2010; Yang & Saslaw 2011; Di 
Porto et al. 2014), and can be therefore used to estimate b by 
comparison with the cosmic variance of the dark matter 

expected from the theory. The galaxy bias computed from the 
cosmic variance, noted bv, encodes the clustering information of 
(Tv and has to be the same that the bias estimated with an in¬ 
dependent method at similar scales. We use the bias computed 
from the two-point correlation function, noted b^, as a bench¬ 
mark because it is a well-tested method widely used in the liter¬ 
ature. The agreement between b,, and b^ implies that the origin 
of the observational cTv is the clustering of galaxies, as desired, 
while the discrepancy reflects deficiencies in our methodology 
to measure c^^,. 

We take advantage of the unique design, depth, and pho¬ 
tometric redshift accuracy of the Advanced, Large, Homoge¬ 
neous Area, Medium-Band Redshift Astronomical (ALHAM¬ 
BRA) survey^ (Moles et al. 2008) to study the clustering sig¬ 
nal encoded in the cosmic variance. We compare the observed 
cosmic variance with the expectations from non-linear theory 
to estimate the galaxy bias bv Then, we compare b^ with the 
b^ reported by Arnalte-Mur et al. (2014, AM 14 hereafter) from 
two-point correlation function analysis using the same data set, 
cosmological parameters, and theoretical expectations. 

The paper is organised as follows. In Sect. 2, we present 
the ALHAMBRA survey and its photometric redshifts and prob¬ 
ability distribution functions, and in Sect. 3 our estimation of 
the cosmic variance and the galaxy bias. We compare our re¬ 
sults with the bias estimated from correlation function analysis 
in Sect. 4. Linally, we present our summary and conclusions 
in Sect. 5. Throughout this paper we use a standard cosmology 
with Dm = 0.27, Da = 0.73, Hq = IQQh km s“^ Mpc^', and 
(Tg = 0.816. To avoid systematics in our results, these cosmo¬ 
logical parameters are the same as in AM14. Magnitudes are 
given in the AB system (Oke & Gunn 1983) and the absolute 
B-band magnitudes are given as Mb - 5 logjg h, even when not 
explicitly indicated. 

2. ALHAMBRA survey 

The ALHAMBRA survey provides a photometric data set 
over 20 contiguous, equal-width (~300A), non-overlapping, 
medium-band optical filters (3500A- 9700A) plus three stan¬ 
dard broad-band near-infrared (NIR) filters (J, H, and K^) over 

' http://www.alhambrasurvey.com 


eight different regions of the northern sky (Moles et al. 2008). 
The survey has the aim of understanding the evolution of galax¬ 
ies throughout cosmic time by sampling a large cosmologi¬ 
cal fraction of the universe, for which reliable spectral energy 
distributions (SEDs) and precise photometric redshifts (zp) are 
needed. The simulations of Benitez et al. (2009), which relate 
the image depth and the accuracy of the photometric redshifts to 
the number of filters, have demonstrated that the filter set chosen 
for ALHAMBRA can achieve a photometric redshift precision 
that is three times better than a classical 4-5 optical broad-band 
filter set. The final survey parameters and scientific goals, as well 
as the technical requirements of the filter set, were described by 
Moles et al. (2008). The survey has collected its data for the 
20 h- 3 optical-NIR filters in the 3.5m telescope at the Calar Alto 
observatory, using the wide-held camera LAICA (Large Area 
Imager for Calar Alto) in the optical and the OMEGA-2000 
camera in the NIR. The full characterisation, description, and 
performance of the ALHAMBRA optical photometric system 
were presented in Aparicio-Villegas et al. (2010). A summary of 
the optical reduction can be found in Cristobal-Hornillos et al. 
(in prep.), the NIR reduction is reported in Cristobal-Hornillos 
et al. (2009). 

The ALHAMBRA survey has observed eight well-separated 
regions of the sky. The wide-held camera LAICA has four 
chips, each with a 15' x 15' held of view (0.22 arcsec pixeU*). 
The separation between chips is 13'. Thus, each LAICA point¬ 
ing provides four distinct areas in the sky, one per chip. Six 
ALHAMBRA regions comprise two LAICA pointings. In these 
cases, the pointings dehne two separate strips in the sky. We 
assumed the four chips in each LAICA pointing to be inde¬ 
pendent subhelds. We summarise the properties of the seven 
helds included in the hrst ALHAMBRA data release in Table 1. 
Currently, ALHAMBRA comprises 48 subhelds of ~ 183.5 
arcmin^, which can be assumed to be independent for cosmic 
variance studies as demonstrated by Lopez-Sanjuan et al. (2014). 

2.1. Bayesian photometric redshifts in ALHAMBRA 

The photometric redshifts used throughout were fully presented 
and tested in Molino et al. (2014), and we summarise their prin¬ 
cipal characteristics below. 

The photometric redshifts of ALHAMBRA were estimated 
with BPZ2.8, a new version of the Bayesian Photometric Red¬ 
shift (BPZ, Benitez 2000) estimator. This is a SED-htting 
method based on a Bayesian inference, where a maximum like¬ 
lihood is weighted by a prior probability. ALHAMBRA re¬ 
lied on the ColorPro software (Coe et al. 2006) to perform 
point spread function (PSL) matched aperture-corrected pho¬ 
tometry, which provided both total magnitudes and isophotal 
colours for the galaxies. In addition, a homogeneous photomet¬ 
ric zero-point recalibration was performed using either spectro¬ 
scopic redshifts (when available) or accurate photometric red¬ 
shifts from emission-line galaxies (Molino et al. 2014). Sources 
were detected in a synthetic F%\AW filter image, noted I in the 
following, defined to resemble the HST/L814W filter. The areas 
of the images affected by bright stars and those with lower expo¬ 
sure times (e.g. the edges of the images) were masked following 
Arnalte-Mur et al. (2014). The total area covered by the current 
ALHAMBRA data after masking is 2.38 deg^ (Table 1). Linally, 
a statistical star/galaxy separation was encoded in the variable 
Stellar_Flag of the ALHAMBRA catalogues, and we kept 
ALHAMBRA sources with Stellar_Flag < 0.5 as galaxies. 

The photometric redshift accuracy, as estimated by compar¬ 
ison with ~ 7200 spectroscopic redshifts (zs), was encoded in 
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Table 1. ALHAMBRA survey fields. 


Field 

name 

Overlapping 

survey 

RA 

(J2000) 

DEC 

(J2000) 

subfields / area 
(no. / deg2) 

ALHAMBRA-2 

DEEP2 (Newman et al. 2013) 

02 28 32.0 

-hOO 47 00 

8/0.377 

ALHAMBRA-3 

SDSS (Aihara et al. 2011) 

09 16 20.0 

-h46 02 20 

8 / 0.404 

ALHAMBRA-4 

COSMOS (Scoville et al. 2007) 

10 00 00.0 

-H02 05 11 

4 / 0.203 

ALHAMBRA-5 

GOODS-N (Giavalisco et al. 2004) 

12 35 00.0 

-h61 57 00 

4/0.216 

ALHAMBRA-6 

AEGIS (Davis et al. 2007) 

14 16 38.0 

-h52 24 50 

8 / 0.400 

ALHAMBRA-7 

ELAIS-Nl (Rowan-Robinson et al. 2004) 

16 12 10.0 

-h54 30 15 

8 / 0.406 

ALHAMBRA-8 

Total 

SDSS (Aihara et al. 2011) 

23 45 50.0 

-h15 35 05 

8/0.375 

48/2.381 


the normalised median absolute deviation (NMAD) of the pho¬ 
tometric vs spectroscopic redshift distribution (Ilbert et al. 2006; 
Brammer et al. 2008), 

C^NMAD - 1-48 X (1) 

where — Zp- Zs and (■) is the median operator. The fraction of 
catastrophic outliers t] is defined as the fraction of galaxies with 
I S, 1/(1 +Zs) > 0.2. In the case of ALHAMBRA, o-nmad = 0.011 
for I < 22.5 galaxies with a fraction of catastrophic outliers of 
rj = 2.1%. We refer to Molino et al. (2014) for a more detailed 
discussion of the ALHAMBRA photometric redshifts. 


The B-band absolute magnitude of a galaxy with observed 
magnitude 1-20 and spectral type T that is located at redshift 
z is denoted as (z, T), which is also provided by BPZ2.8. 
This function condenses the information about the luminosity 
distance, which depends on z, and the k correction, which de¬ 
pends on z and T, needed to translate observed /-band magni¬ 
tudes to absolute B-band magnitudes. Thus, the Mg of the galaxy 
i depends on its z and T, and we have to weight (z, T) with 
the corresponding PDF, (z, T) to obtain , . Formally, we esti¬ 
mate Mg as a function of z as 


Mg^iiz) 


f Mf(z, r)xPDF,(z, r)d7’ 
f PDF,(z, 7’)dr 


+ (/,• - 20 ), 


(5) 


2.2. Probability distribution functions in ALHAMBRA 

This section is devoted to the description of the probability dis¬ 
tribution functions (PDFs) of the ALHAMBRA sources. When 
dealing with photometric redshifts, several studies show that it 
is better to use the full PDF than to use only the best Zp (e.g. 
Femandez-Soto et al. 2002; Cunha et al. 2009; Wittman 2009; 
Myers et al. 2009; Schmidt & Thorman 2013; Carrasco Kind 
& Brunner 2014). The ALHAMBRA PDFs have been success¬ 
fully used to study high-redshift (z > 2) galaxies (Viironen et al. 
2015), to detect galaxy groups and clusters (Ascaso et al. 2015), 
to estimate the merger fraction (Lopez-Sanjuan et al. 2015), or 
to improve the estimation of stellar population parameters (Diaz- 
Garcia et al. 2015). 

The probability that a galaxy i is located at redshift z and has 
a spectral type T is PDF, (z, T). This probability density function 
is the posterior provided by BPZ2.8. The probability that the 
galaxy i is located at redshift z is then 

PDF, (z) = J PDF, (z,/’) dr. (2) 

Moreover, the total probability that the galaxy i is located at zi < 
Z < Z2 is 

r,(zi,Z 2 )= rPDF,(z)dz. (3) 

dzi 

The probability density function PDF, (z, T) is normalised to one 
by definition; in other words, the probability of any galaxy i be¬ 
ing found in the whole parameter space is one. Formally, 


1 = 




PDF,(z)dz = 



PDF,(z, r)drdz. 


(4) 


where /, is the observed /-band magnitude of the source. 

Thanks to the probability functions defined in this section, 
we are able to statistically use the output of current photometric 
redshift codes without losing information. 


2.3. Sample selection 

We focus our analysis on the galaxies in the ALHAMBRA first 
data release^. This catalogue comprises ~ 500k sources and is 
complete (5cr, 3" aperture) for / < 24.5 galaxies (Molino et al. 
2014). 

We perform our study in a given redshift range z e 
Umim^max) with B-band luminosity selected samples. Because 
the Mg of a galaxy depends on its redshift as shown in Sect. 2.2, 
we manage the galaxy samples under study with the function S, 
defined as 


Si(z,Q,Mf) = 


1 , 

0 , 


if Mgj (z) + Qxz< Mf, 
otherwise. 


(6) 


where Mg i (z) is the B-band luminosity of the galaxy i from 
Eq. (5), is the selection magnitude of the sample, and the 
term Q x z accounts for the evolution of the luminosity func¬ 
tion with redshift. This last term ensures that we explore similar 
Lg/L*g luminosities at every z, where L*g is the typical B-band 
luminosity at a given redshift reported by Ilbert et al. (2006). 

Our goal is to compare the galaxy bias derived from the cos¬ 
mic variance with the bias measured from the two-point corre¬ 
lation function analysis performed by AM 14. To minimise any 
possible difference between the two studies, we define the same 
redshift ranges and B-band luminosity selections as AM14. The 
studied samples, defined with Q = 0.6 to mimic AM14 selection, 

^ http://cloud.iaa.es/alhambra/ 
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are summarised in the first two columns of Table 2, and cover a 
B-band luminosity from LbIL*^ - 0.2 to 2 at 0.35 < z < 1.05. 
The LbIL*^ values in the table are those reported by AM14. 

3. Estimation of the gaiaxy bias from the cosmic 
variance 

In this section we present the steps we follow to estimate the 
galaxy bias from the cosmic variance in ALHAMBRA. First, we 
expose the theoretical relation between the cosmic variance and 
b (Sec. 3.1). Then, we detail how to measure the cosmic variance 
from ALHAMBRA data (Sec. 3.2). Finally, we estimate the cos¬ 
mic variance of the dark matter expected from theory (Sec. 3.3). 


3.1. Theoretical link between cr^ and b 


To connect the cosmic variance with the galaxy bias, we follow 
Somerville et al. (2004) and Moster et al. (2011). The mean 
{N} and the variance {N^} - {N}^ in the distribution of galax¬ 
ies are given by the hrst and second moments of the probability 
distribution (V), which describes the probability of counting 
N objects within a volume V. The relative cosmic variance is 
defined as 


2 _ - (Nf 1 


(7) 


The second term represents the correction for the Poisson shot 
noise. The second moment of the object counts is 

(At2) = (Nf + (N}+^-^£ f(|r, - rbl) dV, dV^, (8) 


where ^ is the two-point correlation function of the sample under 
study (Peebles 1980). Combining this with Eq. (7), the relative 
cosmic variance can be written as 

crl^ ^ j^^(\r,-r^,\)dV,dV^. (9) 

We can approximate the galaxy correlation function in Eq. (9) 
using the correlation function for dark matter as ^ — b^ ^dm, 
where b is the galaxy linear bias and we assume that b does not 
depend on scale. With this definition of the correlation function, 
we find that 
a2 r 

^ ^ j fdmdra - TbDdVadVb = (^0) 


where ctj, dm is the cosmic variance of the dark matter. Einally, 
the galaxy bias from the cosmic variance is 



v,dm 


( 11 ) 


We estimate cr,, from ALHAMBRA data in Sect. 3.2, and cr,, dm 
from theory in Sect. 3.3. 



Fig. 1. Number density distribution of Mg' = -18.1 (top panel) and 
Mg' = -20.6 (bottom panel) galaxies at 0.55 < z < 0.85 from the 48 
ALHAMBRA subfields. The star and the black bar mark the median 
and the intrinsic dispersion retrieved by the MLE, respectively. The 
red solid line shows the best MLE solution convolved with the observa¬ 
tional errors and it is independent of the histogram binning. The derived 
relative cosmic variance cr„ is labelled in the panels. A colour version 
of this plot is available in the electronic edition. 


where Nj is the number of galaxies that fulfil the selection, and 
Vj is the cosmic volume probed by the subfield j at Zmin < z < 
Zmax- The index i spans all the galaxies in the subfield; that is, no 
pre-selection of the sources is performed. A similar probabilistic 
approach is used by Viironen et al. (2015) to study the number 
counts of high-redshift galaxies. 

The 48 number densities tij from the ALHAMBRA subfields 
follow a log-normal distribution. 


Bln (n I o") 


n V^( 


exp 


(\nn-p)^ 

2cr^ 


(13) 


where p and cr are the median and the dispersion of a Gaussian 
function in log-space. 


3.2. Measuring the cosmic variance rr^ in ALHAMBRA 


In this section we define the procedure used to measure the cos¬ 
mic variance from ALHAMBRA data. We estimate the num¬ 
ber density of galaxies (n) in the ALHAMBRA subfield j for 
a given B-band luminosity selection and redshift range 
Z ^ [Zmin; ^max) 


Nl 

yj 



PDE,(z)x.S,(z,0.6,Mf)dz, 


( 12 ) 


Bg (n I p, 0 -) 


V^c 


exp 


(n'-p? 

2cH- 


(14) 


where n' - Inn. This log-normal distribution of the num¬ 
ber densities was expected because the distribution of overdense 
structures in the universe is described by a log-normal function 
(e.g. Coles & Jones 1991; de la Torre et al. 2010; Kovac et al. 
2010; Yang & Saslaw 2011). We present two representative ex¬ 
amples in Eig. 1. We checked that the number densities of all the 
samples under study follow a log-normal distribution with an 
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Table 2. Galaxy bias estimated from the cosmic variance in ALHAMBRA. 


Mf - 5 logio h 
(1) 

LbILI 

(2) 

n 

(10-3/r^Mpc-3) 

(3) 

o-v 

(4) 

by 

(5) 

(6) 


0.35 < 

z 

< 0.65 z = 

0.52 CTyjJj^ = 

0.175 


-16.8 

0.16 


32.4 + 0.8 

0.174 + 0.019 

1.00 + 0.11 

1.00 + 0.09 

-17.6 

0.27 


20.7 + 0.6 

0.200 + 0.023 

1.14 + 0.13 

1.16 + 0.11 

-18.1 

0.37 


15.0 + 0.5 

0.216 + 0.025 

1.24 + 0.14 

1.23 + 0.13 

-18.6 

0.50 


10.1 +0.4 

0.228 + 0.027 

1.30 + 0.15 

1.27 + 0.12 

-19.1 

0.73 


6.2 + 0.2 

0.251 +0.031 

1.43 + 0.17 

1.35 + 0.17 

-19.6 

1.06 


3.3 + 0.1 

0.271 +0.036 

1.55 + 0.20 

1.50 + 0.16 

-20.1 

1.52 


1.4+ 0.1 

0.318 + 0.049 

1.81 +0.28 

1.79 + 0.17 

-20.6 

2.40 


0.42 + 0.04 

0.387 + 0.096 

2.21 + 0.55 

2.40 + 0.50 


0.55 < 

z 

< 0.85 z = 

0.74 cr^,_dm = 

0.148 


-17.6 

0.26 


22.4 + 0.6 

0.174 + 0.019 

1.18 + 0.13 

1.22 + 0.17 

-18.1 

0.35 


16.3 + 0.5 

0.194 + 0.022 

1.31+0.15 

1.34 + 0.21 

-18.6 

0.49 


11.0 + 0.4 

0.215 + 0.025 

1.46 + 0.17 

1.52 + 0.24 

-19.1 

0.69 


6.7 + 0.3 

0.250 + 0.029 

1.69 + 0.30 

1.74 + 0.29 

-19.6 

0.97 


3.5+ 0.1 

0.278 + 0.035 

1.88 + 0.24 

2.00 + 0.30 

-20.1 

1.44 


1.5+ 0.1 

0.325 + 0.044 

2.20 + 0.30 

2.30 + 0.40 

-20.6 

2.13 


0.50 + 0.03 

0.422 + 0.072 

2.85 + 0.49 

3.20 + 0.60 


0.75 < 

z 

<1.05 z = 

0.91 CTy,dm ~ 

0.134 


-18.1 

0.34 


17.6 + 0.6 

0.190 + 0.022 

1.41+0.16 

1.37 + 0.15 

-18.6 

0.46 


11.7 + 0.3 

0.198 + 0.023 

1.48 + 0.17 

1.47 + 0.15 

-19.1 

0.66 


7.0 + 0.2 

0.220 + 0.026 

1.64 + 0.19 

1.62 + 0.17 

-19.6 

0.93 


3.5+ 0.1 

0.247 + 0.030 

1.84 + 0.22 

1.90 + 0.30 

-20.1 

1.35 


1.4+ 0.1 

0.330 + 0.044 

2.47 + 0.33 

2.21 + 0.23 

-20.6 

1.95 


0.38 + 0.03 

0.427 + 0.070 

3.18 + 0.52 

2.90 + 0.30 


Notes. Col. (1): Absolute 6-band magnitude at z = 0 used to select the sample with Eq. (6). Col. (2): Median luminosity of the sample in units 
of L* from AM14. Col. (3): Median number density of the sample. Col. (4): Observed cosmic variance. Col. (5): Galaxy bias measured from the 
cosmic variance. Col. (6): Galaxy bias measured from the two-point correlation function by AM14. 


Anderson & Darling (1952) test. This test confirms that the dis¬ 
tribution of the number densities is compatible with a log-normal 
shape and always disfavours a Gaussian distribution. 

There are two origins of the observed cr; the intrinsic disper¬ 
sion (Tint (i.e. the field-to-field variation due to the clustering of 
the galaxies), and the dispersion due to the Poisson shot noise 
CTp. We estimate the Poisson shot noise in each subfield as 

crpj = ^|NJ/VJ (15) 

and apply the maximum likelihood estimator (MLE) presented 
in Lopez-Sanjuan et al. (2014) to measure the median number 
density of the distribution n, the intrinsic dispersion cri„,, and 
their uncertainties. 

Applying Eq. (7) to Pln («), we conclude that the observed 
relative cosmic variance is 

cry=e°'“-l. (16) 

We estimate the cr„ uncertainty by propagating the uncertainty in 
(Tint, as estimated with the MLE (see Lopez-Sanjuan et al. 2014, 
for details). We summarise the median number density n and 
the measured cosmic variance cr^ of the studied samples in Cols. 
3 and 4 of Table 2. We note that the average number densities 
reported by AM14 differ from our values. This discrepancy is 


explained by the number density definition; we estimate a prob¬ 
abilistic number density from Eq. (12) and AM14 assumed the 
galaxies located at their best photometric redshift. 

3.3. Estimation of the dark-matter cosmic variance cr„d„ 

The final ingredient needed to estimate the galaxy bias from 
Eq. (11) is cr^ (jnj, the cosmic variance of the dark matter expected 
from the theory. We calculate CTy^dm in each volume by solving 
the integral in Eq. (9) for dark matter using the code QUICKCV^, 
which is described in Newman & Davis (2002). The code com¬ 
putes the cosmic variance from the dark-matter power spectrum 
using a window function which is 1 inside the interest volume 
and 0 otherwise. We obtain the dark-matter power spectrum at 
each redshift bin using the CAMB software (Lewis et al. 2000), 
including the non-linear corrections of HALOFIT (Smith et al. 
2003). 

The cosmic variance is a particular case of count-in-cell 
statistics, and the procedure used to estimate the galaxy linear 
bias is similar if either spherical or cubic volumes are used in¬ 
stead (see references in Sect. 1). In the next section, we compare 

^ QUICKCV is available at www.phyast.pitt.edu/~janewman/ 
quickcv 
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the galaxy bias by estimated from the cosmic variance with the 
values reported in ALHAMBRA by AM14 from two-point cor¬ 
relation function analysis. They estimate the galaxy bias as 


b( = (17) 

V VVp,dm 

where Wp is the measured projected correlation function and 
vt'p.dm is the projected dark-matter correlation function expected 
from the theory. 


4. Clustering signal encoded in the cosmic variance 

We summarise the studied samples, the measured cosmic vari¬ 
ance cr,„ the inferred galaxy bias from the cosmic variance by, 
and the galaxy bias b^ reported by AM 14 from correlation func¬ 
tion analysis in Table 2. We show the comparison between by 
and b^ in three redshift bins and for several B-band luminosity 
selections in Fig. 2. 

We find that the galaxy bias from the cosmic variance nicely 
agrees with the values from AM 14, with all the estimations con¬ 
sistent at the Icr level. This result has two important implica¬ 
tions. First, it confirms that the cosmic variance measured in AL¬ 
HAMBRA is due to the clustering of galaxies and the method¬ 
ology applied to measure cTy is therefore correct. Thus, the AL¬ 
HAMBRA cTy can be used to characterise the cosmic variance 
affecting pencil-beam surveys, reinforcing the results and the 
parametrisation of (jy for close-pair studies presented by Lopez- 
Sanjuan et al. (2014). Second, we can estimate the galaxy bias 
b from the cosmic variance, complementing the bias measure¬ 
ments from correlation function analysis with an independent 
method. The cr,, is measured in large volumes defined by the ob¬ 
servational strategy of the survey, and it is unaffected by the pho¬ 
tometric redshift precision and peculiar motions as the counts-in- 
cells performed with spherical or cubic volumes (Papageorgiou 
et al. 2012). 

Regarding the observed trends, we find that by increases with 
the B-band luminosity of the sample and, for a fixed selection, 
increases with redshift. These trends have already been found 
in several studies (e.g. Benoist et al. 1996; Norberg et al. 2001; 
Polio et al. 2006; Coil et al. 2006; Coupon et al. 2012; Skibba 
et al. 2014; AM14). The quantitative comparison with the galaxy 
bias measured by other surveys is beyond the scope of this paper. 
Because of the excellent agreement with the AM14 values, we 
refer the reader to AM14 for a detailed discussion of the galaxy 
bias measurements in ALHAMBRA and a comparison with pre¬ 
vious work in the literature. 

In the next sections we explore the impact on the by vs 
b^ comparison of the ALHAMBRA fields explored (Sect. 4.1), 
the power-spectrum model used (Sect. 4.2), the probed scale 
(Sect. 4.3), and the colour of the galaxies (Sect. 4.4). 
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Fig. 2. Galaxy bias b as a function of the B-band luminosity in three 
redshift bins (labelled in the panels). The dots are from this work by 
cosmic variance analysis, b,.. The red crosses are from Arnalte-Mur 
et al. (2014) from the two-point correlation function analysis, b^. The 
dotted line marks unity. A colour version of this plot is available in the 
electronic edition. 


4.1. Dependence on the explored fields 

Arnalte-Mur et al. (2014) conclude that the fields ALHAMBRA- 
4 (COSMOS) and ALHAMBRA-7 (ELAIS-Nl) have signifi¬ 
cantly different clustering properties from the other ALHAM¬ 
BRA fields. Thus, in addition to the global measurements 
presented in the previous section, AM 14 also reported the 
b^ estimated without using the fields ALHAMBRA-4 and 
ALHAMBRA-7. To further check the ALHAMBRA cos¬ 
mic variance, we also estimated the galaxy bias by with¬ 
out the 12 subfields corresponding to the ALHAMBRA-4 and 
ALHAMBRA-7 fields. The results are summarised in Fig. 3. As 


in the global case, the agreement between both estimations of the 
galaxy bias is remarkable and further support that the clustering 
measured by the correlation function is encoded in the estimated 
cosmic variance, as desired. 

4.2. Dependence on the power-spectrum model 

The power spectrum that we used to estimate CTy^iyn accounts 
for the non-linear evolution of the dark matter. The popular 
work of Moster et al. (2011) also uses QUICKCV but with a linear 
power spectrum. We checked that the values from Moster et al. 
(2011) are recovered if linear theory is used, and we estimated 
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the impact of non-linear evolution in the theoretical models. We 
find that the crv^dm from non-linear theory is larger by 15%, 9%, 
and 6% than the cosmic variance from linear theory reported 
by Moster et al. (2011) at z = 0.52, 0.74, and 0.91, respec¬ 
tively. These differences are below the 20% accuracy targeted 
by Moster et al. (2011) and decrease with redshift for two main 
reasons: first, the impact of the non-linear evolution increases 
with the growth of the structures, and therefore decreases with 
z. Second, the non-linear effects are important at scales smaller 
than ~ l/r ' Mpc. The scale probed by each ALHAMBRA sub- 
field increases with z, decreasing the impact of the non-linear 
evolution. 

4.3. Dependence on the probed scale 
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Throughout this paper we assume that the galaxy bias b does 
not depend on scale, but b is a scale-dependent parameter (e.g. 
Heavens et al. 1998; Mann et al. 1998; Cresswell & Percival 
2009). The observational studies find that the bias is nearly con¬ 
stant for scales larger than r ~ 1/;“' Mpc (the linear bias regime), 
with larger b values at smaller scales (e.g. White et al. 2011; 
Parejko et al. 2013). The values of b^ reported by AM14 are 
computed in the linear bias regime, from Mpc to lO/z^* 
Mpc. In the cosmic variance case, we integrate all the scales 
inside the volume of interest, and the measurements of b^ are 
affected by scales smaller than r = l/z ' Mpc. The integration 
in Eq. (9) weights each scale by r^, so the impact of the small 
scales decreases with the area subtended by the ALHAMBRA 
subfields. This subtended area increases from ~ 3.6 x 3.6 
Mpc^ at z = 0.52 to ~ 4.5 x 4.5/i“^ Mpc^ at z = 0.91. We find 
that bv and b^ are similar, implying that b^ is dominated by the 
linear bias regime and that the ALHAMBRA subfields are large 
enough to reduce the scale dependence of the bias. 

We further test the impact of small scales with a toy model. 
In our model, we assumed a scale-dependent bias of the form 


b(r) = 


4-2r, 

2 , 


if r < Ih ' Mpc, 
if r > l/i“' Mpc, 


(18) 


as suggested by observations (White et al. 2011; Parejko et al. 
2013). We applied this bias model to the non-linear power spec¬ 
trum and computed the cosmic variance of this biased power 
spectrum, noted The comparison between and 

o’v.dm provides an effective bias b^^ that should be close to 
bes — 2 if the signal is dominated by the linear bias regime. 
We find fieff = 2.2,2.1, and 2.1 at z = 0.52, 0.74, and 0.91, 
respectively. This confirms that small scales should affect our 
measurements by less than 10%, supporting that b^, is compara¬ 
ble to the linear bias measured with b^ given the geometry of the 
ALHAMBRA survey. 


4.4. Dependence on the colour 

Following the work of AM 14, Hurtado-Gil et al. (2015) study 
the two-point correlation function and the galaxy bias of red and 
blue galaxies with < -18.6 in the ALHAMBRA survey. To 
complete the comparison between b^ and b^, we also estimated 
the galaxy bias of red and blue galaxies with M^®* < -18.6 from 
the cosmic variance. 

To define red and blue galaxies, we take advantage of the 
profuse information encoded in the PDFs. Instead of selecting 
galaxies according to their observed colour or their best spectral 
template, we split each PDF into red templates (T - E/SO), de¬ 
noted PDF'^®'*, and blue templates (T - S/SB), denoted PDF*’*“®. 


4.0 

z 1 1 1 

^ 1 Z 

3.5 

0.55 < z < 0.85 

O by This work 

- 

3.0 

- ™ AM14 

- 

2.5 


A 

2.0 


T 

1.5 

1.0 

. 1 . 

- 


0.0 0.5 1.0 1.5 2.0 2.5 



Lb/L'b 

4.0 

z - \ -1-1-1-\—: 

3.5 

0.75 < z < 1.05 

1 


O by This work 1 

3.0 

- H bj AM14 1 

^ 2.5 

} T 

2.0 

i 

1.5 


1.0 

— 1 ■ 1 ■ 1 ■ 1 ■ ■ 1 ■ 1 ■ 1 ■ ■ 1 ■ 1 ■ 1 ■ ■ 1 ■ 1 ■ 1 ■ ■ ■ ■ 1 ■ 1 ■ ■ ■ ■ 1 ■ 1 ■ 1 - 

_1_1_1_1_1_ 


0.0 0.5 1.0 1.5 2.0 2.5 

Lb/L*^ 


Fig. 3. Galaxy bias fc as a function of the S-band luminosity in 
three redshift bins (labelled in the panels) with the 12 subfields from 
ALHAMBRA-4 and ALHAMBRA-7 removed. The dots are from this 
work by cosmic variance analysis, b,. The red crosses are from Amalte- 
Mur et al. (2014) by two-point correlation function analysis, b^. The 
dotted line marks unity. A colour version of this plot is available in the 
electronic edition. 


Formally, 
PDF(z) = 


PDP®‘‘ (z) H- PDF““® (z) 


'blue . 


/■ 


/ 


PDF(z,E/S0)d7’-H PDF(z,S/SB)d7’. (19) 


In practice, the red templates have T 6 [1,5.5] and the blue 
templates have T e (5.5,11] in the ALHAMBRA catalogues. 
Thus, we can reliably work with colour segregations without any 
pre-selection of the sources (e.g. Lopez-Sanjuan et al. 2015). 

We summarise the basic properties of the samples under 
study, the measured cry and by, and the b^ reported by Hurtado- 
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Table 3. Galaxy bias of galaxies with < -18.6 as a function of colour estimated from the cosmic variance in ALHAMBRA. 


Redshift range 
(1) 

n 

(10-^ Mpc- 
(2) 

"^) o-y 

(3) 

^ v,dm 

(4) 

by 

(5) 

b( 

(6) 

Red templates (E/SO) 

0.35 < z < 0.6 

3.4 + 0.2 

0.349 + 0.046 

0.195 

1.84 + 0.26 

1.69 + 0.08 

0.6 < z < 0.8 

3.0+ 0.2 

0.317 + 0.045 

0.183 

1.78 + 0.26 

1.55 + 0.13 

0.8 <z< 1.0 

2.7 + 0.2 

0.244 + 0.033 

0.163 

1.52 + 0.21 

1.83 + 0.11 

Blue templates (S/SB) 

0.35 < z < 0.6 

7.2+ 0.3 

0.229 + 0.032 

0.195 

1.19 + 0.17 

1.18 + 0.06 

0.6 < z < 0.8 

7.4 + 0.3 

0.224 + 0.031 

0.183 

1.24 + 0.18 

1.13 + 0.07 

0.8 <z< 1.0 

8.4 + 0.3 

0.190 + 0.024 

0.163 

1.18 + 0.15 

1.10 + 0.20 


Notes. Col. (1): Redshift range under study. Col. (2): Median number density of the sample. Col. (3): Observed cosmic variance. Col. (4): 
Dark-matter cosmic variance from the theory. Col. (5): Galaxy bias measured from the cosmic variance. Col. (6): Galaxy bias measured from the 
two-point correlation function by Hurtado-Gil et al. (2015). 



Fig. 4. Galaxy bias b of E/SO (red symbols) and S/SB (blue symbols) 
galaxies with Mg' < -18.6 as a function of redshift. The dots are 
from this work by cosmic variance analysis, The crosses are from 
Hurtado-Gil et al. (2015) by two-point correlation function analysis, 
bf. The dashed lines and the coloured areas mark the error-weighted 
average of b^ and its error, respectively, for red = 1.69 ± 0.14) 
and blue = 1.20 ± 0.10) galaxies. The dotted line marks unity. A 
colour version of this plot is available in the electronic edition. 


Gil et al. (2015) in Table 3. We note that the subfields 
from ALHAMBRA-4 are not included in this analysis because 
Hurtado-Gil et al. (2015) discard this held owing to its peculiar 
clustering properties. As in the general case, agrees with 
at the Icr level both for red and blue galaxies, reinforcing our 
results and conclusions (Fig. 4). We hnd that red galaxies are 
more biased than blue galaxies, with by being compatible with 
a constant bias in the redshift range under study. The error- 
weighted averages at 0.35 < z < 1.0 are b\f‘^ - 1.69 + 0.14 
and = 1.20 ± 0.10. This implies a relative bias of hrei = 
= 1.4 + 0.2, in agreement with previous work that hnd 
byei ~ 1.3 - 1.7 at z < 1 (e.g. Madgwick et al. 2003; Meneux 
et al. 2006; Coil et al. 2008; de la TorTe et al. 2011; Skibba et al. 
2014; Hurtado-Gil et al. 2015). 


5. Summary and conclusions 

We estimate the signihcance of the clustering signal encoded in 
the cosmic variance cTy measured in the ALHAMBRA survey. 
With this aim, we measured the galaxy bias b from the ALHAM¬ 
BRA cosmic variance. The cosmic variance is dehned as the in¬ 
trinsic dispersion (i.e. due to the held-to-held variations) in the 
number density distribution of the 48 ALHAMBRA subhelds. 
The number densities were computed using the full probability 
distribution functions (PDFs) of the ALHAMBRA photometric 
redshifts. We compared the observed cry with the cosmic vari¬ 
ance of the dark matter expected from non-linear theory to esti¬ 
mate the galaxy bias b. 

The galaxy bias from the cosmic variance nicely agrees with 
the bias estimated by AM14 in ALHAMBRA from two-point 
correlation function analysis. This confirms that the cosmic vari¬ 
ance measured in ALHAMBRA is due to the clustering of galax¬ 
ies and the methodology applied to measure cTy is therefore cor¬ 
rect. Thus, the ALHAMBRA cry can be used to characterise the 
cosmic variance affecting pencil-beam surveys, reinforcing the 
results and the parametrisation of the cosmic variance for close- 
pair studies presented by Lopez-Sanjuan et al. (2014). We hnd 
that the bias increases with the Z?-band luminosity of the sample 
and, for a hxed selection, increases with redshift, in agreement 
with previous studies (e.g. Norberg et al. 2001; Coupon et al. 
2012; AM14). Moreover, red galaxies have a larger bias than 
blue galaxies, with a relative bias of hrei = 1.4 + 0.2 (e.g. Madg¬ 
wick et al. 2003; Meneux et al. 2006; Skibba et al. 2014). 

The technique outlined in this paper can be used to esti¬ 
mate the galaxy bias b of high-redshift galaxies and of differ¬ 
ent galaxy populations at z < 1 (Sect. 4.4) simply from the dis¬ 
persion of the observed number densities, complementing the 
bias measurements from correlation function analysis with an 
independent method. These number densities are measured in 
large volumes dehned by the observational strategy of the survey 
rather than in the common spherical or cubic volumes, which are 
affected by the photometric redshift precision and peculiar mo¬ 
tions (Papageorgiou et al. 2012). In addition to its cosmological 
and observational value, the study of the number density distri¬ 
bution can be useful in order to detect systematics and photomet¬ 
ric inhomogeneities in the next generation of large photometric 
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surveys, such as the Dark Energy Survey (DES, Flaugher 2012), 
the Javalambre-Physics of the accelerating universe Astrophysi- 
cal Survey (J-PAS, Benitez et al. 2014), and the Large Synoptic 
Survey Telescope (LSST, Ivezic et al. 2008). 
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